Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

2QQg  2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2008  to  00-00-2008 

4.  TITLE  AND  SUBTITLE 

Dynamics  of  Coronal  Mass  Ejections  (CMEs) 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROTECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory ,4555  Overlook  Avenue 

SW, Washington, DC, 20375 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBIECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

_ _ _  ABSTRACT 

18.  NUMBER  19a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  OS 

unclassified  unclassified  unclassified  Report  (SAR) 

3 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


SIMULATION,  COMPUTING,  AND  MODELING 


Dynamics  of  Coronal  Mass  Ejections  (CMEs) 

J.  Chen,1  P.W.  Schuck,1  and  V.  Kunkel2 
Plasma  Physics  Division 
2 George  Mason  University 

Introduction:  Strong  geomagnetic  storms  cause 
disturbances  in  the  ionosphere  that  can  modify  con¬ 
ditions  in  the  operational  environment  of  the  Navy 
around  the  globe.  If  severe  enough,  they  can  degrade 
or  cause  outages  in  SATCOM  and  GPS  operations, 
impairing  the  Navy’s  ability  to  maintain  environmental 
awareness  and  continuous  connectivity  of  the  widely 
dispersed  fleet.  The  ultimate  cause  of  strong  storms  is 
the  eruption  of  Earth-directed  coronal  mass  ejections 
(CMEs)  at  the  Sun.  CMEs,  which  are  strongly  magne¬ 
tized  plasma  structures,  represent  ejection  of  I016  g  of 
material  at  up  to  2000  km/s  from  the  Sun. 

The  magnetic  geometry  and  the  driving  mecha¬ 
nism  of  CMEs  have  been  two  major  questions  of 
modern  solar  physics.  These  questions  are  closely 
related:  the  driving  forces  acting  on  a  magnetized 
plasma  structure  critically  depend  on  its  magnetic 
geometry  in  three  dimensions.  Unfortunately,  neither 
the  magnetic  field  of  the  eruptive  structure  nor  the 
forces  acting  on  it  can  be  directly  measured.  Under  an 
NRL  basic  research  program  supplemented  by  NASA 
grants,  we  have  developed  a  first-principles  theory  to 
understand  the  3D  magnetic  geometry  and  the  nature 
of  the  force  from  directly  observable  quantities. 

What  is  the  3D  Magnetic  Geometry  of  CMEs? 

With  the  launch  of  NRL’s  Large  Angle  and  Spectro- 
metric  Coronagraph  (LASCO)  instruments  on  the 
Solar  and  Heliospheric  Observatory  (SOHO)  satellite 
in  1995,  CME  acceleration — the  net  accelerating  force 
per  unit  mass — has  been  extensively  documented  in  the 
region  2-30  solar  radii  (Rq).  Our  work  with  LASCO 
data  has  shown  that  the  underlying  magnetic  field  of 
a  CME  is  that  of  a  “flux  rope”:  pictorially  speaking,  a 
bent  coil  of  magnetic  held  lines,  whose  two  ends  are 
connected  to  the  Sun.  Figure  3  shows  a  composite 
image  of  a  CME  (arrow)  observed  by  the  LASCO  C2 
and  C3  telescopes  on  20  August  2000.  In  this  figure,  the 
small  white  circle  at  the  bottom  indicates  the  limb  of 
the  Sun.  The  circular  disk  of  radius  2RQ  corresponds 
to  the  occulter  of  the  C2  telescope.  The  two  legs  of  the 
CME  connected  to  the  Sun  are  separated  by  distance  Sf, 
as  illustrated  in  the  figure. 

Using  a  theoretical  model  previously  developed  at 
NRL,  we  demonstrated  that  the  height  Zmax  above  the 
photosphere  at  which  the  acceleration  of  the  leading 
edge  attains  its  maximum  value  scales  with  the  distance 
Sf  between  the  photospheric  footpoints  of  the  CME’s 
legs  in  such  a  way  that  the  following  inequality  is  satis¬ 
fied:1 


(1) 


This  relationship  was  directly  tested  against  17  well- 
observed  CMEs  and  eruptive  prominences  (EPs), 
where  (i)  the  source  regions  were  visible  so  that  we 
could  determine  the  approximate  distance  (Sf)  between 
the  footpoints  and  (ii)  the  data  contained  sufficient 
numbers  of  height  measurement  to  allow  us  to  deter¬ 
mine  the  maximum -acceleration  height  (Zmax).  Radio 
data  and  data  from  the  Extreme  Ultraviolet  Imaging 
Telescope  (EIT)  on  SOHO  were  used  to  estimate  Sf. 

Each  event  was  analyzed  to  determine  a  pair  of 
numbers,  Zmax  and  Sf.  Figure  4,  adapted  from  Ref.  1, 
shows  the  data  points  for  the  17  events.  The  lines  A  and 
B  are  defined  by  Zmax  =  Sf/ 2  and  Zmax  =1.5  Sf,  respec¬ 
tively,  corresponding  to  the  lower  and  upper  limits  of 
the  inequality  above.  The  plot  shows  that  inequality  (1) 
is  satisfied  by  all  the  events  in  this  sample  consistent 
with  the  uncertainties  in  the  data. 


What  is  the  Force  that  Drives  CMEs?  What  can 
we  deduce  from  the  excellent  agreement  between  the 
theoretical  prediction  and  observational  data?  The 
flux- rope  scaling  law,  inequality  (1),  is  based  on  two 
assumptions:  (i)  the  magnetic  structure  is  that  of  a  flux 
rope  with  footpoints  separated  by  Syand  (ii)  the  struc¬ 
ture  is  driven  by  the  “hoop  force,”  a  particular  form 
of  Lorentz  force  acting  on  curved,  current-carrying 
plasmas,  including  tokamaks  in  laboratories  for  mag¬ 
netic  fusion.  This  force  has  the  unique  property  that 
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where  R  is  the  major  radius  of  the  flux  rope.  The 
right-hand  side  is  maximum  where  curvature  (1 IR)  is 
maximum. 

The  physics  of  this  process  is  schematically  illus¬ 
trated  in  Fig.  5.  Imagine  that  the  lowest  arc  corresponds 
to  the  initial  flux  rope  characterized  by  radius  radius  Rj 
or  curvature  1/Ri  and  footpoint  separation  Sf.  As  the 
flux  rope  rises  with  the  footpoints  held  stationary,  the 
curvature  increases.  Curvature  reaches  its  maximum 
value  when  the  arc  is  a  semi-circle  with  Z  =  R  =  Sfl 2, 
where  Z  is  the  height  of  the  apex.  Thereafter,  the  curva¬ 
ture  monotonically  decreases.  Thus,  if  the  flux  rope  is 
driven  by  the  hoop  force,  the  acceleration,  i.e.,  the  force 
per  unit  mass,  necessarily  peaks  when  the  flux  rope 
reaches  the  semi-circular  position.  The  good  agreement 
shown  in  Fig.  4  implies  that  CMEs  are  consistent  with 
magnetic  flux  ropes  and  expand  under  the  Lorentz 
force.  Physically,  the  result  shows  that  CME  accelera¬ 
tion  is  governed  by  the  footpoint  separation  distance 
Sf. 
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FIGURE  3 

A  LASCO  C2-C3  composite  image  of  a 
CME,  indicated  by  an  arrow,  observed 
on  20  August  2000.  The  white  circle  at 
the  bottom  indicates  the  Sun.  The  bright 
object  on  the  right  is  a  planet,  with  the 
horizontal  rays  caused  by  the  saturated 
instrument. 


Sf(R0) 


FIGURE  4 

Comparison  of  theory  and  data.  Zmax  and  Sf  are 
determined  from  data.  The  point  marked  by  "S" 
denotes  a  numerical  simulation  result.  Data  points 
lying  between  lines  A  and  B  satisfy  the  flux-rope 
scaling  law,  inequality  (1). 


How  Does  a  CME  Propagate?  After  the  initial 
acceleration,  a  CME  continually  expands  away  from 
the  Sun.  The  propagation  through  the  interplanetary 
medium,  however,  could  not  be  observed  until  recently. 
With  the  launch  of  a  pair  of  Solar  Terrestrial  Relations 
Observatory  (STEREO)  satellites  with  the  NRL  suite 
of  instruments  called  Sun-Earth  Connection  Coronal 
and  Heliospheric  Investigation  (SECCHI)2  on  board, 
it  has  now  become  possible  to  quantitatively  study  the 
dynamics  of  CMEs  in  interplanetary  space. 

Figure  6  shows  the  trajectory  of  a  CME  observed 
by  SECCHI  (COR2,  HI  1,  and  HI  2)  and  LASCO  (C2 
and  C3)  in  May  2007.  The  height  of  the  leading  edge 
of  the  CME  measured  from  the  solar  surface  is  shown 
by  various  symbols  corresponding  to  different  instru¬ 


ments.  The  curve  is  the  initial-value  solution  of  the 
theoretical  model.  The  figure  shows  that  the  model  is 
capable  of  producing  a  solution  that  accurately  repro¬ 
duces  the  observed  dynamics  throughout  the  field  of 
view,  from  RQ  to  about  3/4  AU. 

We  have  established  with  strong  quantitative 
evidence  that  CMEs  are  magnetic  flux  ropes  accelerated 
by  the  Lorentz  hoop  force.  Figure  6  represents  the  first 
successful  theoretical  modeling  of  CME  acceleration 
in  the  inner  corona  and  the  subsequent  propagation  in 
interplanetary  space. 

[Sponsored  by  ONR  and  NASA] 


232 


SIMULATION,  COMPUTING,  AND  MODELING 


FIGURE  5 

Flux-rope  scaling  law.  The  arcs  schematically 
represent  a  rising  flux  rope. 
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FIGURE  6 

Theoretical  solution  directly  compared  with  SECCHI  data. 
Measured  height-time  data  are  represented  by  various 
symbols  corresponding  to  different  telescopes.  The  curve  is 
the  "best-fit"  theoretical  solution. 
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